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(57) Abstract: The present invention extends the concept of two dimensional channel estimation to multicarrier multiple input mul- 
tiple output communication systems using multicarrier modulated transmission signals impinging from a plurality of (NT) transmit 
antennas at the receiver side. It is assumed that the transmission signals carry a two dimensional data sequence with embedded pilot 
symbol. As basis of a two-stage channel estimation process the plurality of transmit antennas is divided into disjoint ti 
antenna subsets (Amicrol); In a first stage of channel estimation impinging pilot sequences are separated in relation tc 
antenna subsets (Amicrol) by performing a first stage channel estimation to yield tentative estimates of a channel response ir 
dimension of transmission. In a second stage of channel estimation impinging pilot sequences are separated in relation tc 
subsets (Amicrol); by performing a second stage channel estimation for each antenna in each ti 
subset (Amicrol) to yield an estimation of the channel response. 
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TWO-DIMENSIONAL CHANNEL ESTIMATION FOR MULTICARRIER 
MULTIPLE INPUT OUTPUT COMMUNICATION SYSTEMS 



Field of Invention 

The present invention relates to channel estimation for multiple input multiple 
output communication systems, and in particular to a method and related channel 
estimator for multiple input multiple output communication systems using multi- 
carrier modulation schemes. 

Background Art 

The use of coherent transmission techniques in wireless communication sys- 
tems requires the tracking of mobile radio channels, known as channel estima- 
tion. Since the signals transmitted from multiple transmit antennas are observed as 
mutual interference, channel estimation for multiple input multiple output MIMO 
communication systems is different from the single transmit antenna scenario. For 
multiple input multiple output MIMO communication systems using multi-carrier 
modulation schemes the received signal after multi-carrier demodulation is typi- 
cally correlated in two dimensions, i.e. in time and frequency. In the following, 
orthogonal frequency division multiplexing OFDM will be referred to as one typ- 
ical example for mutli-carrier moduilation schemes. The reason for this is that 
orthogonal frequency division multiplexing OFDM and variants thereof are the 
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most popular multi-carrier modulation schemes. 

Communication systems employing multiple transmit and receive antennas, 
known as multiple input multiple output MEMO communication systems, can be 
used with orthogonal frequency division multiplexing OFDM to improve the com- 
munication capacity and quality of mobile radio systems. For orthogonal fre- 
quency division multiplexing OFDM communication systems with multiple trans- 
mit antennas, such as space-time codes as decibed in A. Naguib, N.Seshadri, and 
A. Calderbank: "Space Time Coding and Signal Processing for High Data Rate 
Wireless Communications", IEEE Signal Processing Magazine, pp. 76-92, May 
2000 or spacial multiplexing, different signals are transmitted form different trans- 
mit antennas simultaneously. Consequently, the received signal is the superposi- 
tion of these signals, which implies challenges for channel estimation. Channel 
parameters are required for diversity combining, if space-time codes are used or 
alternatively for separation of superimposed signals if spatial multiplexing is used. 

Approaches to such channel estimation are described in Y. Li, N. Seshadri, 
and S. Ariyavisitakul: "Channel Estimation for OFDM Systems with Transmitte 
Diversity in Mobile Wireless Channels", IEEE Journal of Selected Areas on Com- 
munications, Vol 17., pp. 461-470, March 1999 and Y. Li: "Simplified Channel 
estimation for OFDM Systems with Multiple Transmit Antennas", IEEE Transac- 
tions on Wireless Communications, Vol 1., pp. 67-75, January 2002, in particular a 
channel estimation scheme for orthogonal frequency division multiplexing OFDM 
with multiple transmit antennas based on the difcrete fourier transform DFT. 

Further, the estimators based on the least squares LS and minimum mean squared 
error MMSE criterion for OFDM-MIMO systems have been systematically de- 
rived in Y. Gong and K. Letaief: "Low Rank Channel Estimation for Space-Time 
Coded Wideband OFDM Systems", Proc. IEEE Vehicular Technology Conference 
(VTC2001-F&11), Atlantic City, USA, pp. 722-776, 2001. Related solutions deal 
with one dimensional approaches where a known pilot OFDM symbol is followed 
by L data bearing OFDM symbols. This scheme is applicable for a quasi-static 
environment where the channel does not change significantly during L OFDM 
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symbols, i.e. indoor systems such as wireless local area networks WLAN. 

To accomodate some mobility, the receiver may switch to decision directed 
channel estimation during the reception of the L data bearing OFDM symbols, as 
suggested in Y. Li, N. Seshadri, and S. Aviyavisitakul: "Channel Estimation for 
OFDM Systems with Transmitte Diversity in Mobile Wireless Channels", IEEE 
Journal of Selected Areas on Communications, Vol 17., pp. 461-470, March 
1999 and Y. Li: "Simplified Channel estimation for OFDM Systems with Mul- 
tiple Transmit Antennas'", IEEE Transactions on Wireless Communications, Vol 
1., pp. 67-75, January 2002. However, decision directed channel estimation which 
uses prior decisions of data symbols as pilot symbols, is significantly more com- 
plex than channel estimation schemes relying on pilots only. 

Further, for OFDM-based systems with one transmit antenna, two dimensional 
channel estimation utilizing a scattered pilot grid can be employed, which satisfy 
the sampling theorem in time and frequency. For pilot-symbol aided channel esti- 
mation PACE known pilot symbols , are multiplexed into the data stream. Interpo- 
lation is used to obtain the channel estimate for the information carrying symbols. 
PACE for single carrier systems was introduced in J.K. Cavers: "An analysis of Pi- 
lot Symbol assited Modulation for Rayleigh Fading Channels", IEEE Transactions 
on Vehicular Technology, Vol. VT-40, pp. 686-693, November 1991. In R. Nils- 
son, O. Edfors, M. Sandell, and P. Boerjesson: "An Analysis of Two-Dimensional 
Pilot-Symbol Assisted Modulation for OFDM", Proc. IEEE Intern. Conf. on 
Personal Wireless Communications (ICPWC'97), Mumbai (Bombay), India, pp. 
71-74, 1997 and P. Hoeher, S. Kaiser, and P. Robertson: "Two-Dimensional Pilot- 
Symbol-Aided Channel Estimation by Wiener Filtering", Proc. IEEE Intern. Conf. 
on Acoustics, Speech, and Signal Processing.(ICASSP'97), Munich, Germany, pp. 
1845-1848, 1997 two-dimensional 2D filtering algorithms have been proposed for 
pilot-symbol aided channel estimation PACE. However, such a 2D estimator struc- 
ture is generally too complex for practical implementation. 

To reduce the complexity, separating the use of time and frequency correlation 
has been proposed in P. Hoeher, S. Kaiser, and P. Robertson: "Pilot-Symbol-Aided 
Channel Estimation in Time and Frequency", in Proc. Communication Theory 
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Mini-Conference (CTMC) in conjunction with IEEE Global Telecommunications 
Conference (GLOBECOM'97), Phoenix, USA, pp. 90-96, 1997. This combined 
scheme, termed double one-dimensional (2 x ID) pilot-symbol aided channel esti- 
mation PACE, uses separate Wiener filters, one in frequency direction and one in 
time direction. 

Another approach to reduce the computational complexity is based on a trans- 
formation which concentrates the channel power to a few transform coefficients. 
Estimators based on the discrete Fourier transform DFT have the advantage that a 
computationally efficient transform in form of the FFT does exist, and that DFT- 
based interpolation is simple. In Y. Li: "Pilot-Symbol- Aided Channel Estimation 
for OFDM in Wireless Systems", IEEE Transactions on Vehicular Technology, 
Vol. 49, pp. 1207-1215, July 2000, the approach based on discrete Fourier trans- 
form DFT-based pilot-symbol aided channel estimation PACE was extended to 
two dimensional pilot-symbol aided channel estimation PACE by using the two 
dimensional FFT for the single antenna case. 

However, a major problem for extending the approach to two dimensional chan- 
nel estimation utilizing a scattered pilot grid to multiple input mulitple output 
MIMO communication systems is that the limitation of the number of transmit 
antennas which can be separated by a certain number of pilots. The minimum 
number of pilot symbols Np which are required to estimate N T channel impulse 
responses (CIR) each of which having Q taps has been shown in Y Gong and K. 
Letaief: "Low Rank Channel Estimation for Space-Time Coded Wideband OFDM 
Systems", Proc. IEEE Vehicular Technology Conference (VTC '2001 -Fall), At- 
lantic City, USA, pp. 722-776, 2001 to be 

N^>N T Q (1) 
However, this means that the number of pilots required for channel estimation 
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grows with the number of transmit antennas Nr. 

Summary of Invention 

In view of the above, the object of the present invention is to extend the concept 
ot two dimensional channel estimation to MBVIO systems. 

According to the present invention this object is achieved through a method of 
two dimensional channel estimation for multiple input multiple output transmis- 
sion systems using multicarrier modulated transmission signals impinging from a 
plurality of transmit antennas and carrying a two dimensional data sequence with 
embedded pilot symbols. In a first step the plurality of transmit antennas is di- 
vided into disjoint transmission antenna subsets. In a second step impinging pilot 
sequences are seperated in relation to transmission antenna subsets by performing 
a first stage channel estimation to yield tentative estimates of a channel response 
in a first dimension of transmission. In a third step impinging pilot sequences are 
seperated in relation to antennas in transmission antenna subsets by performing 
a second stage channel estimation for each antenna in each transmission antenna 
subset to yield an estimation of the channel response. 

An important advantage of the present invention is increased flexibilty in chan- 
nel estimation. The reason for this is that by dividing the separation task of the 
superimposed transmission signals in time and frequency direction, a more effi- 
cient usage of the pilot symbols is possible. Hence, either the number of required 
pilot symbols may be reduced or the performance can be improved. 

In view of the above, another important advantage of the present invention is the 
increase in the number of transmit antennas which can be estimated with a certain 
number of pilot symbols through application of a two stage channel estimaiton 
approach. 

Yet another important advantage of the present invention is that the two stage 
channel estimation approach allows for tracking of channel variations even at high 
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Doppler frequencies. This is a prerequisite to support high velocities of mobile 
users and therefore to enable truely mobile multiple input multiple output MIMO 
communication systems. 

According to a preferred embodiment of the present invention the first stage 
channel estimation is performed using pilot sequences arranged as a two dimen- 
sional grid of pilot symbols, wherein pilot symbols used for first stage channel esti- 
mation depend on the first dimension of transmission only and pilot symbols used 
for second stage channel estimation depend on a second dimension of transmis- 
sion only. Preferably, pilot sequences are expressed in a product form for achiev- 
ing seperability of pilot sequences in the first dimension of transmission and the 
second dimension of transmission. 

An advantage of this preferred embodiment of the present invention is that uti- 
lization of a scattered pilot grid allows for efficient use of pilot symbols. Further, 
by matching the pilot spacing in time and frequency to the worst case channel char- 
acteristics higher mobile velocities can be supported with respect to conventional 
one dimensional schemes. 

According to another preferred embodiment of the present invention, for a pilot 
spacing having a value of one the first stage channel estimation and/or the second 
stage channel estimation is achieved in a non-interpolating manner through yield 
of tentative estimates in relation to pilot symbol grid positions in the dimension of 
estimation. Alternatively, for a pilot spacing having a value larger than one the first 
stage channel estimation and/or the second stage channel estimation is achieved in 
an interpolating manner through yield of tentative estimates for all data sequence 
grid positions in the dimension of estimation. 

An important advantage of this preferred embodiment is flexible support of 
different two dimensional pilot grids. In other words, the present invention may 
be flexibly applied using any type of pilot spacing, both, in frequncy and time 
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direction the application of suitable interpolation techniques. 

Further preferred embodiments of the present invention relate selection of first 
dimension of transmission for the first channel estimation stage and the second 
channel estimation stage - i.e., in frequence direction or in time direction - and 
further to the selection of channel estimation domain— i.e., frequency domain 
channel estimation or time domain channel estimation. Here, according to the 
present invention any combination of dimension of transmission and channel esti- 
mation domain is supported. 

The free selectability of dimension of transmission and channel estimation do- 
main is further reason for flexibility of the channel estimation approach according 
to the present invention. It enables optimal consideration of multi-carrier related 
transmission parameters, selected pilot grid structure, and also application of com- 
putationanally most suitable channel estimation techniques. 

According to another preferred embodiment of the present invention the chan- 
nel estimation approach is applied to a celluar communication system with a fre- 
quency reuse factor of one such that base stations and related anntenna arrays form 
the plurality of transmit antennas and such that transmission antenna subsets and 
related transmission antennas are defined in relation to this plurality of transmit 
antennas. 

An important advantage of this preferred embodiment of the present invention 
is the application of the two stage channel estimation techiques as outlined above 
to distributed antennas. In particular, it allows to handle a situation where a mobile 
user roams at a cell border. While data bearing symbols can be protected against 
interference using a channel code or spreading, this is not possible for pilot sym- 
bols. According to the present invention through appropriate definition of subset 
in relation to cells in the celluar communication system. 

According to yet another preferred embodiment of the present invention there is 
provided a computer program product directly loadable into the internal memory 
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of a channel estimator for estimating multiple input multiple output transmission 
channels in two dimensions comprising software code portions for performing the 
steps of the method of two dimensional channel estimation according to the present 
invention when the product is run on a processor of the channel estimator. 

Therefore, the present invention is also provided to achieve an implementation 
of the inventive method steps on computer or processor systems. In conclusion, 
such implementation leads to the provision of computer program products for use 
with a computer system or more specifically a processor comprised, e.g., in a chan- 
nel estimator for estimating multiple input multiple output transmission channels 
in two dimensions. 

The programs defining the function of the present invention can be delivered 
to a computer/processor in many forms, including, but not limited to information 
permanently stored on non-writeable storage media, e.g., read only memory de- 
vices such as ROM or CD ROM discs readable by processors or computer I/O 
attachments; information stored on writable storage media, i.e. floppy discs and 
hard drives; or information convey to a computer/processor through communica- 
tion media such as local area network and/or telephone networks and/or Internet 
or other interface devices. It should be understood, that such media when carry- 
ing processor readable instructions implementing the inventive concept represent 
alternate embodiments of the present invention. 



Description of Drawing 

In the following the best mode and preferred embodiments of the present in- 
vention will be explained with reference to the drawing in which: 

Fig. 1 shows a schematic diagram of an OFDM based multiple input multiple 
output MIMO communication system for explanation of the system model under- 
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lying the present invention; 

Fig. 2 shows a schematic diagram illustrating OFDM modulation and demodu- 
lation, respectively; 

Fig. 3 shows a scattered pilot grid suitable for two dimensional channel estima- 
tion according to the present invention; 

Fig. 4 shows a schematic diagram of a channel estimator for estimating mul- 
tiple input multiple output transmission channels of mutlicarrier communication 
systems according to the present invention; 

Fig. 5 shows a flowchart of operation of the channel estimator shown in Fig. 5; 

Fig. 6 shows a schematic diagram illustrating the principle of 2xlD channel 
estimation underlying the present invention; 

Fig. 7 shows a schematic diagram of a channel estimator for estimating mul- 
tiple input multiple output transmission channels of mutlicarrier communication 
systems according to the present invention, wherein channel estimation is per- 
formed in frequency direction first; 

Fig. 8 shows a further scattered pilot grid suitable for two dimensional channel 
estimation according to the present invention; 

Fig. 9 shows a further scattered pilot grid corresponding to an digital video 
broadcast DVB-T application and suitable for two dimensional channel estimation 
according to the present invention; 

Fig. 10 shows a schematic diagram of a channel estimator for estimating mul- 
tiple input multiple output transmission channels of mutlicarrier communication 
systems according to the present invention, wherein channel estimation is per- 
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formed in time direction first; 

Fig. 11 shows a schematic diagram of an estimator stage adapted to achieve 
channel estimation in the time domain according to the present invention; 

Fig. 12 shows a schematic diagram of a further estimator stage adapted to 
achieve channel estimation in the time domain according to the present invention; 
and 

Fig. 13 shows an application of the two stage channel estimation approach 
according to the present invention to a cellular communication system with a fre- 
quency reuse factor of one. 

Description of Best Mode and Preferred Embodiments 

In the following, the best mode and preferred embodiments of the present in- 
vention will be explained with reference to the drawing. Initially, some basic con- 
siderations underlying differential multiple-length transmit diversity and related 
diversity reception will be explained for a better understanding of the present in- 
vention. 

Basic Considerations 

The present invention addresses pilot-symbol aided channel estimation PACE 
schemes for multi-carrier multiple input multiple output MIMO communication 
systems which are based on the insertion of a pilot grid. Typically for a multi- 
carrier multiple input multiple output MIMO communication system, pilot-symbol 
aided channel estimation PACE - also referred to as PACE in the following - is 
applied across subcarriers in frequency direction and over several multi-carrier 
transmission symbols - e.g., OFDM symbols - in two dimensions, resulting in two 
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dimensional 2D PACE. 

The present invention as explained in the following is not restricted to a par- 
ticual type of multi-carrier multiple input multiple output MEMO communica- 
tion system, and may be applied, e.g., to orthogonal frequency divsion multiplex- 
ing OFDM, discrete multitone transmission DMT, filtered multitone transmission 
FMT, or biorthogonal frequency division multiplexing BFDM. 

Another multi-carrier multiple input multiple output MIMO communication 
system is multi-carrier code divsion mutliple access MC-CDMA where spread- 
ing in frequency and/or time direction is introduced in addition to the multi-carrier 
modulation. Yet another multi-carrier multiple input multiple output MIMO com- 
munication system is a multi-carrier code divsion mutliple access MC-CDMA sys- 
tem with a variable spreading factor, namely variable spreading factor orthogonal 
frequency and code division multiple access VSF-OFCDM. 

From the above, it- should be understood that pilot symbol aided channel es- 
timation PACE may be applied to all multi-carrier multiple input multiple output 
MIMO communication systems operating with transmission signals being corre- 
lated in two dimensions. Therefore, all these multi-carrier multiple input multiple 
output MIMO communication systems may employ the different embodiments of 
the present invention as explained in the following. 

System Model 

Fig. 1 shows a schematic diagram of an OFDM based multiple input multiple 
•output MIMO communication system - also referred to as OFDM system in the 
following - for explanation of the system model underlying the present invention. 
Further, Fig. 2 shows a schematic diagram illustrating OFDM modulation and de- 
modulation in correspondence to the OFDM based multiple input multiple output 
MIMO communication system shown in Fig. 1. 

As shown in Fig. 1, for the considered OFDM-based MIMO system, one 
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OFDM modulator is employed on each transmit antenna, while OFDM demod- 
ulation is performed independently for each receive antenna. The signal stream 
is divided into N c parallel substreams, typically for any multi-carrier modulation 
scheme. The i til substream, also referred to as subcarrier in the following, of the 
I th symbol block, named OFDM symbol in the following, is denoted by X tti . An 
inverse DFT with Nppr points is performed on each block, and subsequently the 
guard interval having N GJ samples is inserted to obtain x ijU . After D/A conversion, 
the signal x(t) is transmitted over a mobile radio channel with response h(t, r). 

As shown in Fig. 1, the received signal at receive antenna v consists of super- 
imposed signals from N T transmit antennas. Assuming perfect synchronization, 
the received signal of the equivalent baseband system impinging at receive antenna 
v at sampling instants t — [n + lN sym }T S pi is in the form 



= y M ([n + W sym ]T spl ) = / h^\t,r)-x^(t-r)dr + n(t) 



where x w {t) denotes transmitted signal of transmit antenna /x after OFDM 
modulation, n(t) represents additive white Gaussion noise, and N sym — Nfft + Nqi 
accounts for the number of samples per OFDM symbol. 

As shown in Fig. 2, at the receiver the guard interval is removed and the infor- 
mation is recovered by performing an DFT on the received block of signal samples, 
to obtain the output of the OFDM demodulation Y e>i . The received signal at receive 
antenna v after OFDM demodulation given by 



where X w t>i and denotes the transmitted information symbol and the 

channel transfer function (CTF) of transmit antenna p, at subcarrier % of the I th OFD. 
symbol, respectively. The term N it i accounts for additive white Gaussian noise 
AWGN with zero mean and variance N 0 . It is assumed that the transmitted signal 
consists of L OFDM symbols, each having N c subcarriers. 




t=[n+£N s 



(2) 
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Channel Characteristics 

While in the following, a way to model channel characteristics will be ex- 
plained, it is important to note that the present invention is not restricted in any 
way through such a model. To the contrary the present invention is applicable to 
actually existing mobile radio channels. For the application of the present inven- 
tion it is only of relevance that mobile radio channels are band limited. Preferably, 
they should also be be limited in time, which practically is the case in view of 
maximum delay of mobile radio channel. Further, preferably different transmit 
antennas and receive antennas should be uncorrelated. 

The present invention considers a time- variant frequency selective fading chan- 
nel. The number of non-zero taps is typically smaller or equal to the maximum 
delay of the channel, Q 0 < Q. The channel impulse response CIR between trans- 
mit antenna n and receive antenna v is defined by 

Qo 

h^\t,T) = J2h^m-t(r-T^) (3) 

9=1 

where v) (t) and u) are the complex amplitude and delay of the q th channel 
tap. 

According to the present invention it is assumed that the Q 0 channel taps and all 
antennas are mutually uncorrelated. The channel taps h^' u) (t) are zero-mean com- 
plex independent identically distributed (i.i.d.) Gaussian random variables. Due to 
the motion of the vehicle h l £' v) (t) will be time-variant caused by the Doppler ef- 
fect. The g th channel tap hf' v) (t) is a wide sense stationary WSS Gaussian process, 
being band-limited by the maximum Doppler frequency u max . 

Further, it is commonly assumed that the channel impulse response CIR is 
approximately constant during one OFDM symbol, so the time dependency of 
the CIR within one OFDM symbol can be dropped, i.e. hf tV) (t) « for 
t 6 [£T, (i+l)T\. Although this is strictly true only for time-invariant channels, 
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this assumption seems to be most often justified in practice, and a good system 
design should ensure that the OFDM symbol duration is sufficiently short. 

Further, the channel transfer function CTF of equation (2), is the the Fourier 
transform of the CIR hP ,v) {t, t). Sampling the result at time t = lT sym and fre- 
quency f = i/T, the CTF at subcarrier % of OFDM symbol I becomes 

Qo 

Hft* = H ^(£T sym ,i/T) = J2KT e ~ j2 * Tqi/T (4) 

q-l 

where T sym = (iNfcr + N GI )T spl and T = iV FFT T spi represents the OFDM symbol 
duration with and without the guard interval, respectively. The matrix form of 
equation (4) is given by 

H'/"'"' = T^hf ^ (5) 

where T 1 "'^ represents the transformation matrix which transforms h' e i,1,v) into 
the frequency domain, defined by 

{T^ v) } g . = exp (-j2ir ^ ; 0 < i < N c - 1 , 1 < q < Qo (6) 

of dimension jV c x Q 0 . 

Further, if the guard interval is longer than the maximum delay of the channel, 
i.e. N G i > Q, where Q > Q 0 denotes the total number of channel taps, the 
orthogonality at the receiver after OFDM demodulation is maintained, and the 
received signal of equation (2) is obtained. 

Assuming the fading at the receiver antennas is mutually uncorrected, the chan- 
nel estimation according to the present invention will be performed independently 
for each antenna. Therefore, in the description to follow the marker for the receive 
antenna v is omitted. 

Further, according to the present invention a channel is defined to be sample 
spaced if the tap delays r q are multiples of the sample instant T spl , i.e. 

r q = T spl P g ; l<q<Qo, P q G {0, 1, • ■ • , Q} (7) 
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where fi g is an arbitrary integer, equal to or larger than zero. In this case H^- v) 
in equation (4) becomes the DFT of the ClRh l [^\ 

In view of the above, Hf- U) can be expressed in matrix notation 

H?"' 10 = Fhf ' v) (8) 
where F represents the DFT-matrix of dimension N G i x N c , defined by 

{F} n . = e -P*"i/*for . 0<i<N c -l, 0<n<N G i-l (9) 

While for a real channel the tap delays r q will not be multiples of the sample 
duration T sp i, however, for many applications the sample spaced channel model 
does approximate a real channel sufficiently well. 

Principle of Pilot Symbol Aided Channel Estimation for OFDM-based MIMO 
Systems 

Fig. 3 shows a scattered pilot grid suitable for two dimensional channel estima- 
tion according to the present invention. 

As shown in Fig. 3, pilot aided channel estimation PACE is based on periodi- 
cally inserting known symbols, termed pilot symbols in the data sequence. If the 
spacing of the pilot symbols is sufficiently close to satisfy the sampling theorem, 
channel estimation and interpolation for the entire data sequence is possible. 

For two dimensional pilot aided channel estimation in the sense of the present 
invention it must be taken into account that for OFDM the fading fluctuations are 
in two dimensions, in time and frequency. In order to satisfy the two-dimensional 
sampling theorem, the pilot symbols are therefore scattered throuout the time- 
frequency grid, which yields a two-dimensional pilot grid. Another reason for 
selecting scattered pilot grids is to maximize spectral efficiency. 
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Description of Pilot Grids in Two Dimensions 



To formally describe a regular grid in the 2D plane according to the present 
invention the total number of pilots transmitted in one frame is defined to N P = 
Np ■ Ng, with A£ = N c /D f and Np = L/D t being the number of pilots in fre- 
quency and time direction respectively. Here, the following notation is used: given 
a matrix describing a 2D structure X, the subsets which describe the dimension 
corresponding to the frequency and time directions are denoted by X' and X", 
respectively. 

Denoting the pilot symbol of subcarrier % and OFDM symbol I by p = [i, £} T , 
any regular 2D grid for use in the framework of the present invention is be de- 
scribed by 



where p = [z, £} T denotes the index of the I th and I th pilot in frequency and time 
direction, respectively, and g 0 defines a pilot grid offset. 

For multiple input multiple output MIMO communication systems - also re- 
ferred to as MIMO system in the follwoing - every transmission signal at transmit 
antenna may to use its own pilot grid. This enables the receiver to separate the 
superimposed transmission signals from different transmission antennas. 

To describe pilot symbol-assisted channel estimation in the sense of the present 
invention there is defined a subset of the received signal sequence containing only 
the pilots, {Xffi = {X^}, with Gp = [i, if € Q, sampled at a D f times lower 
rate z = [i/D f \ in frequency direction, and at a D t times lower rate I = [£/D t \ 
in time direction, respectively. As a general convention, variables describing pilot 
symbols will be marked with a ~ in the following description. 

Further, without restricting scope of protection, it may be assumed that the 




(10) 
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pilots X w i:i are chosen from a PSK constellation, so \Xfl\ = 1. 

For the particular example shown in Fig. 3 the structure of the pilot grid is 
defined by G which is 

r - \ Df 0 " 

L 0 ^. 

where D f and D t denotes the so-called pilot spacing in frequency and time, 
respectively. For the pilot grid shown in Fig. 3 the particular values are D f = 5 
and D c = 5. 

It should be noted that before transmission, the pilots {Xf2} ma Y be multiplied 
by an outer pilot sequence {X 0[l } which is identical for all transmit antennas to 
yield the transmitted pilot sequence 



Further, the outer pilot sequence {X 0/ J is chosen to have a low peak to power 
average ratio in the time domain and/or to have good correlation properties for 
synchronization. 

As shown in Fig. 2, at the receiver the cyclic prefix is removed and an FFT 
is performed to yield the received signal after OFDM demodulation. Assuming 
perfect synchronization, the received signal Y t>i of equation (2) is obtained. For 
channel estimation the received signal at the pilot positions are demultiplexed from 
the data stream, and after removing the outer pilot sequence, by dividing through 
Xo-„ the received pilot is obtained according to 

N T 

where Gp is defined in (10). 
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Basic Consideration for 2x1 D-Pilot Assisted Channel Estimation for Multiple 
Input Multiple Ouput Communication Systems 

Fig. 4 shows a schematic diagram of a channel estimator for estimating mul- 
tiple input multiple output transmission channels of mutlicarrier communication 
systems according to the present invention. 

As shown in Fig. 4, the channel estimatior 10 comprises a first estimator stage 
12 and a second estimator stage 14. Further, the channel estimator comprises a 
transmission antenna subset memory 16. 

Fig. 5 shows a flowchart of operation of the channel estimator shown in Fig. 4. 

As shown in Fig. 5, operation of the channel estimator according to the presen 
invention relies on a method where in a firsat step S10 the plurality of transmit 
antennas into disjoint transmission antenna subsets. Then, in a step S 12 impinging 
pilot sequences in relation to transmission antenna subsets are seperated by per- 
forming a first stage channel estimation to yield tentative estimates of a channel 
response in a first dimension of transmission. Then, in a step S14 impinging pilot 
sequences in relation to antennas in transmission antenna subsets by performing 
a second stage channel estimation for each antenna in each transmission antenna 
subset to yield an estimation of the channel response. 

Fig. 6 shows a schematic diagram illustrating the principle of channel estima- 
tiong according to the present invention. 

As shown in Fig. 6 and outlined above, in the step S12 channel estimation 
is performed in one dimension, yielding tentative estimates for all subcarriers of 
these OFDM symbols. These tentative estimates are then used in step S14 as new 
pilots, in order to estimate the channel for the entire frame. 

As shown in Fig. 6, the second stage does not only interpolate between OFDM 
symbols having pilots, but it does also improve the accuracy of the tentative esti- 
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mates. 

Further, according to the present invention Either channel estimation in fre- 
quency direction or time direction may be performed first. Reference to the case 
where channel estimation in frequency direction is performed first will be made by 
2 x ID-PACE type I in the following. Further, reference to the case where channel 
estimation in time direction is performed first will be mde by 2 x ID-PACE type II 
in the following. 

To formally describe the problem, the received pilot of OFDM symbol IDt is 
considered at the (z\D/) th subcarrier 

Y loau s = E X ^ H ^o f + N io^ f t = {l,2,---,L/D t } 

i = {1,2, ••• ,N C /D f } 

where X Ul) i Dt - iDj and H^\ D denotes the transmitted pilot symbol and the 
channel transfer function (CTF) of transmit antenna fi, at subcarrier % = iDf of the 
£ = lDt h OFDM symbol, respectively. 

Further, it is assumed that the CTF varies in the i and in the % variable, i.e. in 
time and in frequency. The term Nt Dt ,iDj accounts for additive white Gaussian 
noise AWGN. L represents the number of OFDM symbols per frame, N c is the 
number of subcarriers per OFDM symbol, D / and Dt denote the pilot spacing in 
frequency and time, and N T is the number of transmit antennas. 

The overall object underlying the present invention is to estimate H { [J for all 
{I, i, /i} within the frame. 

To achieve this object, generally for MIMO channel estimation, in addition to 
channel estimation and interpolation it is also necessary to separate the impinging 
signals from Nt transmit antennas. Dividing the signals corresponding to the N T 
transmit antennas into N TX subsets, yields Nn groups each of which having N T 2 
signals, such that N T = N T \N T2 . In other words, for a MIMO system having iV T 
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transmit antennas, according to the present invention signals corresponding to N T2 
transmit antennas into one set, to form N T1 subsets. 

In other words, the basic concept underlying the present invention is to divide 
the separation task into two stages, in the way that in step S12 channel estimation 
is performed in one dimension, at OFDM symbols I — £D t , yielding tentative 
estimates for all subcarriers of these OFDM symbols. The second step S14 uses 
these tentative estimates as new pilots, in order to estimate the channel for the 
entire frame. 

As shown in Fig. 6, the second stage estimation does not only interpolate be- 
tween OFDM symbols having pilots, but it does also improve the accuracy of the 
tentative estimates. Therefore, the different embodiments of the present invention 
as explained in detail in the following have significantly reduced complexity while 
there is little degradation in performance. 

According to Fig. 6 channel estimation in frequency direction is performed 
first, to reverse the order such that channel estimation in time direction is per- 
formed first is straightforward. In the following reference will be made to the case 
where channel estimation in frequency direction is performed first by 2 x ID-PACE 
type I. to the contrary, reference to the case where channel estimation in time di- 
rection is performed first will be made to as 2 x ID-PACE type II. 

Considering channel estimation in time direction, it is proposed to use the pilots 
\Yi t i, Y Dt>i , • ■• , Y DtN n+i^ e Q, in order to estimate HfJ. Further, for smoothing 
type filtering it is proposed to use past as well as future pilots to estimate Hfj, this 
means 1 < I < D t Np. Therefore, for smoothing an estimate cannot be obtained 
until all pilots have been received, which requires buffering of A£ = D t N'p — I 
OFDM symbols. 

As will be shown in more detail in the follwoing, the alternative is to use pre- 
diction type filtering where I > D t Np. In this case only past pilots are used for 
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channel estimation in time direction. 

Obviously, prediction type filtering does not require any buffering, however, 
the performance with respect to smoothing degrades. For channel estimation in 
frequency direction, on the other hand, all pilots of one OFDM symbol are be- 
ing received together, so no buffering is required. However, the accuracy of the 
channel estimates typically degrades near the band edges. 

2 x ID-PACE type I 

Fig. 7 shows a schematic diagram of a channel estimator for estimating mul- 
tiple input multiple output transmission channels of mutlicarrier communication 
systems according to the present invention, wherein channel estimation is per- 
formed in frequency direction first. 

According to the present invention it is prposed to extend 2 x ID-PACE to 
OFDM-based MIMO channel estimation. For MIMO channel estimation it is 
necessary to separate the impinging signals from N T transmit antennas. Let the 
MIMO system having N T transmit antennas be denoted by set the A Further, 
according to the present invention the set of N T transmit antennas is devided into 
N Tl subsets An C A, with ^ = {1, • • • , N T1 }. Each subset contains N T2 trans- 
mit antennas, such that N T = N T1 iV T 2. 

In other words, for a MIMO system having N T transmit antennas, according 
to the present invention it is proposed to group the signals corresponding to N T2 
transmit antennas into one set, to form N T i subsets of A. Without restricting 
scope of protection, one may assume that all sets have the same number of transmit 
antennas, and the subsets A^ are disjoint, i.e. each transmit antenna can only be 
in one set. 

The concept underlying me present invention is to divide the separation task 
into two stages, in the way that we first separate a subset of the N T \ < N T signals 
together with channel estimation in the first estimation stage. The resulting signal 
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Zf t ° will be a superposition of iV r2 = N T /N T i signals. 



In the second estimation stage the remaining N T2 superimposed signals are 
seperated for each of the N TX signals of the first estimation stage, together with 
channel estimation in the second dimension, to yield the estimate of the frequency 



Fig. 7 illustrates the basic idea for type I of the proposed scheme. 

It should be noted, that the buffer shown in Fig. 7 is used in order to apply 
smoothing type filtering in time direction. However and as outlined above, the 
present invention is applicable to, both, smoothing and prediction type filtering, 
so the buffer shown in Fig. 7 is optional and depends on the particular channel 
estimation algorithm being used. 

Fig. 8 and 9 show pilot sequence designs which may be used to support the 
two stage approach for OFDM-based MEMO channel estimation accordance to the 
present invention. Here, the pilot grid shown in Fig. 9 corresponds to the DVB-T 
pilot grid according to ETSI EN 300 744, V 1.4.1 (2001-01). Further standards - 
however, not to be considered as restricting scope of protection - would be IEEE 
802.1 la or ETSI TS 101 475 HIPERLAN/2. 

In a more genral sense, in order to extend the pilot sequence design for an 
OFDM-based system with multiple transmit antennas, according to a preferred 
embodiment of the present invention the pilot sequence of transmit antenna fj, is 
defined by {X^}. 

It is proposed to choose a pilot sequences that can be expressed in the product 
form 



response H^. 



Xfi = X^> ■ 



fjt = Ai 2 + Nt2 ■ Oi - 1) , Mi = {1> • • • > n ti) 
ju 2 = {!,••■ ,iV T2 } 



1 < i < N£ , 1 <£ < Np 



(12) 
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where X{^ ] and are the pilot symbols for the first and second stage, re- 
spectively. 

It should be noted that the pilot symbol of the first stage X^ 0 only depends 
on the subcarrier index i, while the pilot symbol of the second stage X^ only 
depends on OFDM symbol I 

This notation implies that {X^ 0 } is identical for all OFDM symbols, indepen- 
tent of I. Correspondingly, the sequence {X^} which accounts for the pilots of 
one subcarrier, is also independent of the subcarrier index i. 

Preferably, the pilot sequences {X ( £ l) } and {X ( 2- 2) } are chosen from orthogonal 
designs, e.g., Walsh sequences or phase shifted sequences. 

Substituting the proposed 2D pilot sequence of (12) into the received pilots in 
(11), the following is obtained 

Y u = X^ l) Xtf Hfl + Nt2 ■ (M1 " 1)] + % 

N Tl N T2 

= X^^2 frg + ^-^'^+Nii (13) 

where is the frequency response of transmit antenna \i. 
Further, the received pilot sequence of subset is given by 

= ^Xg* H l g + r,Ta ' llll - 1)) (14) 

In view of the above, the task of the first estimation stage is to estimate Zf^; 
that is to separate the Nti groups, and then to estimate and interpolate the channel 
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in frequency direction. 

It should be noted that the pilot sequence is constant with respect to the 
subcarrier index i. This means that Zf? is a superposition of N T 2 waveforms ftp 
multiplied with a constant phase term X%*\ 

Further, the outputs of the first estimation stage are subsequently used as inputs 
for the second estimation stage shown 7. 

As shown in Fig. 7, in the second estimation stage the channel is estimated 
in time direction to separate the remaining N T2 signals per subset to yield the 
estimate of the frequency response H { [\. 

Using this framework, according to the present invention it is iiroposed to apply 
available one dimensional schemes for OFDM-based MIMO systems to perform 
channel estimation. 

2 x ID-PACE type II 

Fig. 10 shows a schematic diagram of a channel estimator for estimating mul- 
tiple input multiple output transmission channels of multicarrier communication 
systems according to the present invention, wherein channel estimation is per- 
formed in time direction first. 

As shown in Fig. 10, the major difference of 2 x ID-PACE type II over 2 x ID- 
PACE type I is that that the separation iV ri subsets in the first estimation stage is 
performed in conjunction with channel estimation in time direction. This yields 
for the pilot symbols of 2 x ID-PACE type II: 

*3 = • > A* = to + N T 2 ■ fa - 1) , m = {1, • • • , N Tl } 

1 < i < #p , 1 < I < Nj! 

(15) 



24 



WO 2004/073276 



PCT/EP2003/001495 



From the above, it may be seen that the pilot sequences are very similar to 
2 x ID-PACE type I in equation (12), the only difference is that the subcarrier 
index i and the OFDM symbol index I are exchanged. Substituting the proposed 
2D pilot sequence of equation (15) into the received pilots in equation (11), the 
received pilot sequence can be represented according to 

N Ti N T2 

y u = j2 % + Nt2 ' (Mi " 1)} +% ( 16 ) 

Hi=l /t 2 =l 

where again the subcarrier index i and the OFDM symbol index £ are exchanged 
with respect to (13). 

This allows the separation of the N Ti subsets in time direction. The received 
pilot sequence of subset is defined by 

Zff = J2^T Hf{ + NT2 ' ( " 1 - 1)) (17) 

In view of he above and as shown in Fig. 10, the task of the first estiamtion 
stage is to estimate Z^; that is to separate the N T i groups, and then to estimate 
and interpolate the channel in time direction, i.e. over the I variable. 

As shown in Fig. 10, if smoothing type filtering is used, A£D t OFDM symbols 
need to be buffered in order to estimate Z^ 1 '. For the second estimation stage it 
is proposed estimate the channel in frequency direction to separate the remaining 
Nt2 signals per subset. 

Comparison between 2 x ID-PACE type I and type II 
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It can be shown that the performance of both schemes type I and type II is iden- 
tical for most implementations. However, the computational complexity, in terms 
of number of multiplications required for channel estimation may differ. The ac- 
tual computational complexity very much depends on the system parameters, the 
pilot grid structure and the channel estimation algorithm which is used. Further- 
more, the DSP or hardware architecture may favor one scheme. 

More importantly however, does the selection of a certain pilot grid rule out the 
implementation of a particular scheme. For the pilot grid shown in Fig. 8 and Fig. 
9 both schemes type I and type II can be applied. However, this is not generally 
the case. In order to motivate this problem consider the following grid 



shown in Fig. 8. For such a grid structure it is more appropriate to employ 2 x ID- 
PACE type I, since the pilots in frequency direction are placed along a line. On 
the other hand, successive pilots in time direction are shifted one subcarrier apart, 
which makes it difficult to employ 2 x ID-PACE type II. 

Accordingly, choosing the grid shown in Fig. 9, which is described by 

Gdvb = , . 

1 4 

to employ 2 x ID-PACE type II is more appropriate. The grid according to Fig. 
9 has been chosen for the DVB-T standard, so for DVB-T channel estimation ac- 
cording to a preferred embodiment of the present invention in time direction should 
be performed first if 2 x ID-PACE is to be used. 

Applications of 2 x ID-PACE type I 

To describe pilot symbol-assisted channel estimation according to the present 
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invention it is proposed to define a subset of the received signal sequence contain- 
ing only the pilots, {Y~} = {Y$}, with {i, £} G Q, sampled at a £>/ times lower 
rate % — [i/Df\ in frequency direction, and at a D t times lower rate I — [£/D t \ 
in time direction, respectively. 

In the following two channel estimation techniques are described for 2 x ID- 
PACE type I, the first is to estimate the channel in the frequency domain by Wiener 
filter interpolation. The second approach is to transfer the received signal into the 
time domain. 

Frequency Domain Channel Estimation 

First Stage Wiener Filtering 

For 2 x ID-PACE type I, the first step is to estimate the channel in the frequency 
direction. In vector notation, the received pilot sequence of OFDM symbol i be- 
comes 

N T i 

Y' t = ^Xx^zf' + N', GC^* 1 (18) 

where Z^ (Ml) represents the received pilot sequence of subset A^ whose entries 
are defined in equation (14). 

Further, the transmitted pilot sequence, the received pilot sequence of subset 
A^, and the additive noise term, of OFDM symbol I transmitted from antenna \x x , 
are given by 

X^ 1 ' = diag(A^---,X^) G C^ x ^ 
27 
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Further, according to the present invention channel estimation in the frequency 
domain is preferably performed with an FIR interpolation filter, which can be ex- 
pressed for the first stage in frequency direction 

N' p N' p N T i 

i=l i=l Mi=l 

= W ,(pl> |i]Y£ (19) 



It should be noted that in general, the filter W? (All) [i] depends on the location of 
the desired symbol, i.e. the subcarrier index i. This means that not only for every 
transmit antenna but also for every subcarrier a different filter is required. 

The optimum approach to estimate is to use a Wiener interpolation filter 
for W' CMl) [i]. A Wiener filter minimizes the means squard error MSE between the 
pilots sequence and the desired response. It is also known as the minimum MSE 
or equivalently MMSE estimator. 

Further, in order to generate a Wiener filter knowledge about the channel statis- 
tics are required, which are described by the covariance matrix. The covariance 
matrix of the pilots in frequency direction is defined by = E\Y'- t Y 1 ^]. The 
entry of the m th row and n th column of the covariance matrix is given by 

{ R WL,„ = E [%-n %- m] = E [ Y ^-D f n >?<- D , m ] {I, i} € G(20) 

Further, define the cross correlation vector between the frequency response of 
the desired sample Zff and the pilots Y£. The m th entry of the cross correlation 
vector R'%£[i] - Effiff Y'- } can be expressed as 

{KM m = E = E [z?fY;- {i __ m)Df ] ; i = |i/£>,J(21) 

The quantities according to equation (20) and equation (21) are necessary to 
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evaluate the Wiener interpolation filter. The optimum solution in the MMSE sense 
may be determined using the Wiener-Hopf according to 

W^f^R^.R^- 1 (22) 



Second Stage Wiener Filtering 

In view of the above, performing equation (19) for the iV T i subsets and for each 
sucarrier is processed further in the second stage. Here, filtering and interpolation 
in time direction yields the frequency response estimate which is achieved in the 
form 

N'J, 

= X>f M ^ %' > I* = & + N T2 • (/il - 1) = {1, • • • , N T } 

i=\ 

(23) 

where wf^l^At] represents the FIR interpolation filter of the second stage of. 
OFDM symbol I with filter delay Al = D t At. 

A positive Al imposes a time delay of Ai symbols at the receiver output. Then 
the estimation filter is a smoothing type filter. On the other hand, setting Al = 0 
specifies a linear prediction receiver without an induced time delay due to channel 
estimation, at the expense of a somewhat poorer estimate of the CIR. 

It should be noted that the best performance is generally achieved if A£ = 
Np/2, i.e. the symbol to be estimated is in the middle of the pilot sequence which 
is used for estimation of that symbol. Therefore, for A£ = N'J,/2 there are Np/2 
future and past pilot symbols involved. 

Further, in matrix notation equation (23) becomes 

% = w" w [t,M] Zf* v = to + Nn-fa-l) = {!,■■. ,N T } (24) 
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where w" w [£, A£] = A£], • ■ • , w'^, ifl) [£, A£]] is the channel estima- 

tion filter of stage two, and ZJ /tMl) = [Zff, ■ ■ , Z l p.] T denotes the block Np 
outputs of stage one of subcarrier i. 

In analogy to the first estimation stage, the optimum approach to estimate Hf\ 
is to use a Wiener interpolation filter for w" w [£, A£]. The covariance matrix of the 
first stage outputs in time direction is defined by R"-- 0 = E[ Z"^ 1 ' Z- /0n) ]. The 
entry of the m th row and n th column of the covariance matrix is given by 

{R'g ) } mjn = E [zfjT. %> /] (25) 

Further, define the cross correlation vector between the frequency response of 
the desired sample HfJ and the outputs of the first stage Z- /<Ml) . The m th entry of 
the cross correlation vector Bf'^[e, A£] = E[H^ Mti zf M)H ] can be expressed 
as 

{R"^[£,Ai]} m = E-[h% %^/] ; £= [£/D t ] +A£ (26) 

The quantities in equation (20) and equation (21) are necessary to evaluate the 
Wiener interpolation filter. The optimum solution in the MMSE sense is derived 
using the Wiener-Hopf equation according to 

W M [£, A£] = R"^[£, A£] ■ R"ig 1,_1 (27) 

Time Domain Channel Estimation 

Fig. 1 1 shows a schematic diagram of. an estimator stage adapted to achieve 
channel estimation in the time domain according to the present invention. 

As shwon in Fig. 11, an alternative approach to determine the frequency re- 
sponse Hi t i according to the present invention is to estimate the channel impulse re- 
sponse (CIR), h' iin w in the time domain first. In terms of the CIR vector of OFDM 
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symbol I impinging from the ^ th transmit antenna, h' e W = [h' 1 ^, ■ ■ • , h'T&Y-- the 
received pilot sequence of subset becomes 

= £ X%>> F h' t W V = li 2 + N T2 . - 1) (28) 

where F denotes an ivp-point DFT matrix defined by 

{F} n ,i = e~ j2 * ni M ; 0 < i < J\£ - 1 , 0<n<Q-l (29) 

It is assumed that the CIR is time limited to Q < N GI samples. Strictly speak- 
ing this is only true for a sample spaced channel model explained above. For a 
non-sample spaced channel model where the DFT of H e>i is not time limited, so 
oversampling is required in order to avoid aliasing. In this case, Q accounts for the 
number of significant taps. 

Transforming Z^" 1 ' into the time domain by an iVp-point IDFT yields 

ir = F*zf 0 = F* £ F hf v = »2 + N T 2 ■ (/* - 1) 

Ai2=l 

N T2 

= E*r h ~f (30) 

JU2 = 1 

This implies that two stage channel estimation for OFDM-based MIMO sys- 
tems can be employed in the time and frequency domains in the same way. 

The received pilot sequency after OFDM demodulation is given by 

Y' t = Xi Z' e + N't = Xx F Nti z' e + N' t (31) 
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where 



= I z / > • • • , z ! 



f,---,f) e c NtiN 'p xNtiQ 



In order to estimate z^ (Ml) according to the present invention the received pilot 
sequency after OFDM demodulation is transformed into the time domain. 

Further, for time domain channel estimation we choose to pre-multiply Y[ by 
the transmitted pilot sequence X\ and then to transform the result into the time 
domain via an Np— point IDFT, that is 

= (x : F,v ri ) K Y^ = Dft; e C*"« xl 

= D^D^ + D^N; (32) 
where the definition = X x F^ T1 has been introduced. 

In the following two basic estimator structures to determine z' (Ml) according to 
the present invention will be described, namely the least squares LS estimator and 
the minimum mean squared error MMSE estimator. 



Least Squares LS Estimator 



Provided the the inverse of D'f D^- does exists, the least squares (LS) estimator 
may be determined according to 

z'ls, = (D'fD'i)" 1 ^ (33) 
= (D'f ~ l D'f Y'i 
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Since the estimator depends on the transmitted signal, the pilot sequence should 
be properly chosen. The LS estimator exists if D £ - is full rank, unfortunately this 
is not always the case. A necessary condition for the LS estimator to exist is 

N^>N T1 Q (34) 



According to a preferred embodiment of the present invention, two times over- 
sampling provides a good trade-off between minimizing the system overhead due 
to pilots and optimizing the performance, i.e. Np » 2N Ti Q. It is assumed that 
Ngi > Q, i.e. the guard interval is longer than the maximum delay of the channel. 

Further, it should be noted that the LS estimator for more than one transmit 
antenna does only exist in the time domain. 

Minimum Mean Squared Error MMSE Estimator 

The MMSE estimator is given by an FIR filter which is for time domain channel 
estimation 

N' P N' p Nti 

m=l m=l /n = l 

= w /(fll, [n] £' t n = {1, • • • , N c } 

(35) 

The Wiener filter is determined by the Wiener-Hopf equation 

w^V] = R'ft'N • R^T 1 n = {1, • • • , N c ] (36) 



In general, the Wiener filter w /(,n) [n], depends on the location of the desired 
symbol n. In order to generate the MMSE estimator knowledge of the correlation 
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matrices and R'^fra] are required 

R^ = E^i'f] = D'f %D;- eC QN ™* QNn (37) 
' = D'f X'iBfa X'f DJ + N 0 D'f DJ 
= D'f D J R' g2 D'f + iV 0 D'f 

and 

R'S'N = eC 1 ^ 1 (38) 

= R'^[n]D'fD;- 

Further, the covariance matrix of z^- is denoted by R~ = E{z'^ z'^ H } with 
dimension QN TX x QN T \. Furtermore, R'^'fn] is row n+ - 1)Q of R~. The 
covariance matrix in the time domain R~ 5 is related to the covariance matrix in the 
frequency domain R^g by 

R» = F^R'^F^ (39) 

Here, it is assumed that the fading of different transmit antennas is uncorrelated. 
It should be noted that while the LS estimator requires D^ to be full rank, while 
the MMSE estimator requires invertability of R^ as seen from equation (36). For 
this to hold, however, D'^ need not to be full rank. Thus, the MMSE estimator can 
exist even if Np < N Ti Q. 

For the case that D£ is full rank, the inverse of D'f D~ does exist. Then the 
Wiener filter of equation (36) can be simplied to 

w ,,w) [n] = R'Sfln]- ^ g +(D'fD',-) _1 iVo] (D'f Dj)" 1 

2g> = R'^in] ■ |Ri 2 + (D'f Dj) N Q j * ■ ^> (40) 

Further it should be noted that according to the present invention the separation 



34 



WO 2004/073276 



PCT/EP2003/001495 



of the N T1 signals, which is performed by the LS estimator, can be separated from 
the filtering task. 

As outlined above, the MMSE estimator is in general dependent on the choice 
of the pilot symbols. However, choosing orthogonal pilot sequences X^ 1 ' the 
estimator becomes independent on the transmitted pilots. For orthogonal pilots 
X^ Xj" 0 = I • 5 Mjm , where 5 MiJ71 denotes the Kronecker symbol, it can be shown 
that D'f DJ = jVfJ. 

Therefore, the above LS estimator in equation (33) as well as the MMSE es- 
timator in equation (40) can be grossly simplified, since the matrix inversion re- 
quired in equation (33) and equation (40) become straightforward. 

Further, it can be seen from equation (33) and equation (40) that the estimator 
has become independent of the chosen pilot sequence, which does significantly 
simplify the filter generation. 

Fig. 1 1 shows a block diagram of channel estimation and interpolation in the 
time domain using orthogonal pilot sequences. First, the received pilot sequence 
is split into N Ti branches and each branch pre-multiplied by XY" 1 '. Then each 
branch of the received pilot sequence is transformed to the time domain. By means 
of - optional - filtering and/or windowing the channel is estimated in the time 
domain. Zero padding of the first stage estimate extends its lenght to N c samples. 
The estimate of the CTF of an entire OFDM symbol (pilots and data), is obtained 
by an iV c -point FFT of the CIR estimate 

Tl = F^ T1 z' or Z' 1 " 1 ' = F z' (Ml) (41) 

where ¥ Nti is a N Tl N c x N T iN c block diagonal matrix, consisting of N Ti blocks 
of Appoint DFT matrices F. 

As shown in Fig. 1 1, the output of the first stage Z /<Ml> can be fed to the second 
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stage estimator in equation (23). 

Alternatively, z' (Ml) may be fed into (23) to yield the CIR estimate hfj which 
is then transformed into the frequency domain with N Tl FFTs. Moreover, channel 
estimation of the second stage: in the Doppler domain is possible, that is Z'" 11 ' or 
z /(Ml) are transformed into the Doppler domain using equivalent algorithms as in 
the time domain. 

Fig. 12 shows a schematic diagram of an estimator stage adapted to achieve 
channel estimation in the time domain according to the present invention. 

From Fig. 12 it can be seen that the channel estimation of the second stage 
may also be performed with DFT-interpolation cooresponding to the first estima- 
tion stage discribed above. However, it may be computationally more efficient to 
perform the second estimation stage in the time domain as well, i.e., before zero 
padding and the iV c -pomt EFT. 

As shown in Fig. 12, the filtering of the second stage itself remains unaffected 
as the correlation function in frequency and time - i.e. the first and second dimen- 
sion of transmission - are mutually independent. Therefore, of Wiener filtering is 
chosen for the second stagem equations (23) and (27) referred to above may still 
be used, the only difference being that the input becomes instead of 
and the output is hfj^ instead of H\^, where n is the sample index in the time 
domain. After the second stage filtering each of the overall N T outputs - i.e., N T2 
outputs per subset - are transformed back into the frequency domain. 

Application to a Cellular System with a Frequency Reuse Factor of One 

Fig. 13 shows an application of the two stage channnel estimation approach 
according to the present invention to a cellular communication system with a fre- 
quency reuse factor of one. 

As shown in Fig. 13, instead of having an antenna array with Nt antenna 
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elements, the proposed scheme can be applied to distributed antennas as well. E.g., 
an application is to employ 2xlD-PACE to a celluar system with a frequency reuse 
factor of one. In a scenario where the mobile user is at the cell border, the user 
will receive the desired signal from one base station and one or several interfering 
signals from other base stations. 

Further, one may assume that each base station has N T2 antenna elements. 
While the data bearing symbols can be protected against interference using a chan- 
nel code or by spreading, the pilot symbols cannot be protected in this way. Ac- 
curate channel estimation, however, is most important for the system to work effi- 
ciently. One solution is to boost the pilots; this however will increase the interfer- 
ence to users served by other base stations, and thus limits the system capcity. 

According to the present invention, 2 x ID-PACE can be applied to this scenario 
as follows: the base stations form jV^i subsets, each subset having an antenna array 
with N T2 antenna elements, to form an resulting array of N T — NtiNt2 elements. 
This would require inter-cell synchronization. 

Abbreviations 



AWGN Additive white Gaussian noise 
CIR Channel impulse response 
CTF Channel transfer function 
DFT Discrete Fourier transform 
FFT Fast Fourier transform 
GI Guard interval 

IDFT Inverse discrete Fourier transform 
IFFT Inverse fast Fourier transform 
LS Least squares 
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MIMO Multiple input multiple output, generally a system having several transmit 
and receive antennas 

MMSE Minimum mean squared error 

MSE Mean squared error 

OFDM Orthogonal frequency division multiplexing 
PACE Pilot-symbol aided channel estimation 



List of Commonly Used System Parameters 



N FFT EFT length. 
N c Number of subcarriers. 
N G : Number of samples of the guard interval. 
L Number of OFDM symbols per frame. 
T OFDM symbol duration. 
T sp i Sample interval, given by T spi = T/Nfft- 

T sym Total OFDM symbol duration including the guard interval T spl - T+N G iT sp i. 

Qo Number of non-zero channel taps. 
Q Total number of channel taps. 

N R Number of receive antennas. 

N T Number of transmit antennas. 
Nti Number of subsets. 

Nt2 Number of transmit antennas within one subset. 
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Np Number of pilots in frequency direction. 
Np" Number of pilots in time direction. 

Np Number of total pilots used for channel estimation (Np = NpNp). 
Df Pilot spacing in frequency. 
D t Pilot spacing in time. 

List of Commonly Used Variables 

Transmitted OFDM symbol of the p th transmit antenna of OFDM symbol t 
at subcarrier %. 

x^ii) Transmitted signal of the ju th transmit antenna after OFDM modulation. 

Y$ Received OFDM symbol of the v th receive antenna of OFDM symbol £ at 
subcarrier i. 

y w (t) Received signal of the u th receive antenna at time t before OFDM demodu- 
lation. 

Hfl v) Channel transfer function (CTF) of the u th receive antenna arriving from the 
/i th transmit antenna of OFDM symbol i at subcarrier i. 

h^' v) (t) Channel impulse response (CIR) of the u th receive antenna arriving from the 
(j} h transmit antenna. 

hfn } Sampled CIR of the V th receive antenna arriving from the /i* transmit an- 
tenna, at the n th sample of OFDM symbol I. 

l\i Sample of AWGN with zero mean and variance N 0 of OFDM symbol I at 
subcarrier i. 

n(t) Realization of a AWGN process at time t before OFDM demodulation. 
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Claims 



1 . Method of two dimensional channel estimation for multiple input mul- 
tiple output transmission systems using multicarrier modulated transmission sig- 
nals impinging from a plurality of (N T ) transmit antennas and carrying a two di- 
mensional data sequence with embedded pilot symbols, comprising the steps: 

dividing the plurality of transmit antennas into disjoint transmission antenna 
subsets (A^); 

separating impinging pilot sequences in relation to transmission antenna subsets 
(Afr) by performing a first stage channel estimation to yield tentative estimates of 
a channel response in a first dimension of transmission; 

seperating impinging pilot sequences in relation to antennas in transmission 
antenna subsets C4 Ml ) by performing a second stage channel estimation for each 
antenna in each transmission antenna subset C4 M1 ) to yield an estimation of the 
channel response. 

2. . Method according to claim 1, characterized in that the first stage 
channel estimation is performed using pilot sequences arranged as a two dimen- 
sional grid of pilot symbols, wherein pilot symbols used for first stage channel 
estimation depend on the first dimension of transmission only and pilot symbols 
used for second stage channel estimation depend on a second dimension of trans- 
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mission only. 

3. Method according to claim 2, characterized in that pilot sequences 
are expressed in a product form for achieving seperability of pilot sequences in the 
first dimension of transmission and the second dimension of transmission. 

4. Method according to claim 3, characterized in that pilot sequences 
are selected from orthogonal designs. 

5. Method according to claim 4, characterized in that orthogonal de- 
signs are selected from a group comprising Walsh sequences and phase shifted 
sequences. 

6. Method according to one of the claims 1 to 5, characterized in that 
the first dimension is frequency dimension. 

7. Method according to claim 6, characterized in that pilot sequences 
are expressed in the product form according to 

= X™ ■ X™ , 11 = 1.12 + N T2 • (Aii - 1) 

= {!,-■■ ,N Tl } 
\x 2 = {l,--- ,N T2 } 

wherein 

I referes to a vector of data symbols in the two dimensional data sequence ar- 
ranged in time direction; 



41 



WO 2004/073276 



PCT/EP2003/001495 



% is a sub-carrier index in frequency direction; 

X^S* represent pilot symbols for the first stage channel estimation in the first di- 
mension of transmission; 

X { 2- 2) represent pilot symbols for the second stage channel estimation in a sec- 
ond dimension of transmission; and 

Nn is the number of transmission antenna subsets CA Ml ), Nt2 is the number of 
antennas in the transmission antenna subsets (A^), and a product of the num- 
ber of transmission antenna subsets and the number of antennas in the transmis- 
sion antenna subsets equals the number iV T of the plurality of transmit antennas 
N T = N Ti ■ N T2 . 

8. Method according to one of the claims 6 or 7, characterized in that for 
a pilot spacing having a value of one first stage channel estimation and/or second 
stage channel estimation is achieved in a non-interpolating manner through yield 
of tentative estimates in relation to pilot symbol grid positions in the dimension of 
estimation. 

9. Method according to one of the claims 6 or 7, characterized in that for 
a pilot spacing having a value larger than one first stage channel estimation and/or 
second stage channel estimation is achieved in an interpolating manner through 
yield of tentative estimates for all data sequence grid positions in the dimension of 
estimation. 

10. Method according to one of the claims 7 to 9, characterized in that 
the first stage channel estimation is achieved using the product form of pilot sym- 
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bols and a representation of received pilot symbols according to 



, + N T2 ■ (mi - 1)> _)-jV"|- 




wherein 

^2=1 

is a representation of a received pilot sequence to be estimated in relation to the 
transmission antenna subset ; 

H { p. is the frequency response of transmit antenna /r, and 

Nij is noise superimposed during transmission. 

1 1 . Method according to one of the claims 6 to 10, characterized in that 

the first stage channel estimation is achieved in frequency domain. 

12. Method according to one of the claims 6 to 10, characterized in that 

the first stage channel estimation is achieved in time domain. 

13. Method according to claim 1 1 or 12, characterized in that the sec- 
ond stage channel estimation is achieved in frequency domain. 

14. Method according to claim 11 or 12, characterized in that the sec- 
ond stage channel estimation is achieved in time domain. 

15. Method according to one of the claims 1 to 5, characterized in that 
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the first dimension is time dimension. 

1 6. Method according to claim 1 5, characterized in that pilot sequences 
are expressed in the product form according to 

X$> = X[f-X%\ A4 = At2 + iV r2 -( Ml -l), mi = {!,••• ,N T i} 

/i 2 = {!,■■■ ,N T2 } 

wherein 

[ referes to a vector of data symbols in the two dimensional data sequence ar- 
raged in time direction; 

i is a sub-carrier index in freqeuncy direction; 

X[' l . l) represent pilot symbols for the first stage channel estimation in the first di- 
mension of transmission; 

X^ 2> represent pilot symbols for the second stage channel estimation in a sec- 
ond dimension of transmission; and 

Nti is the number of transmission antenna subsets (y4 Ml ), Nt2 is the number of 
antennas in the transmission antenna subsets CA W ), and a product of the num- 
ber of transmission antenna subsets and the number of antennas in the transmis- 
sion antenna subsets equals the number Nt of the plurality of transmit antennas 
N T = N T i • N T2 . 

17. Method according to one of the claims 15 or 16, characterized in 
that for a pilot spacing having a value of one first stage channel estimation and/or 
second stage channel estimation is achieved in a non-interpolating manner through 
yield of tentative estimates in relation to pilot symbol grid positions in the dimen- 
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sion of estimation. 



18. Method according to one of the claims 15 or 16, characterized in 
that for a pilot spacing having a value larger than one first stage channel estima- 
tion and/or second stage channel estimation is achieved in an interpolating manner 
through yield of tentative estimates for all data sequence grid positions in the di- 
mension of estimation. 



19. Method according to claim one of the claims 16 to 18, characterized 
in that the first stage channel estimation is achieved using the product form of 
pilot symbols and a representation of received pilot symbols according to 

Ml = l M2 = l 

wherein 

N T2 

is a representation of a received pilot sequence to be estimated in relation to trans- 
mission antenna subset ; 

Hf-. is the frequency response of transmit antenna /j; and 
Nfj is noise superimposed during transmission. 



20. Method according to one of the claims 15 to 19, characterized in 
that the first stage channel estimation is achieved in frequency domain. 



21. Method according to one of the claims 15 to 19, characterized in 
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that the first stage channel estimation is achieved in time domain. 

22. Method according to claim 20 or 21 , characterized in that the sec- 
ond stage channel estimation is achieved in frequency domain. 

23. Method according to claim 20 or 21 , characterized in that the sec- 
ond stage channel estimation is achieved in time domain. 

24. Method according to claim 1 1 or 20, characterized in that the first 
stage channel estimation in frequency domain is achieved through finite impulse 
response filtering of received pilot symbols to yield tentative estimates (Z 1 ^) of 
received pilot sequences (Zff) in relation to transmission antenna subsets (A^). 

25. Method according to claim 13 or 22, characterized in that the sec- 
ond stage channel estimation in frequency domain is achieved through finite im- 
pulse response filtering of tentative estimates (Zff) in relation to antennas in trans- 
mission antenna subsets (A H ) to yield estimates of the channel reponse (Hffi. 

26. Method according to claim 24 or 25, characterized in that the sec- 
ond stage channel estimation in frequency domain is executed through finite im- 
pulse response filtering with induced time delay for smoothing type filtering. 

27. Method according to claim 24 or 25, characterized in that the sec- 
ond stage channel estimation in frequency domain is executed through finite im- 
pulse response filtering without induced time delay for predictor type filtering. 

28. Method according to one of the claims 24 to 27, characterized in 
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that the finite impulse response filtering is a Wiener type filtering. 

29. Method according to claim 12 or 21, characterized in that the first 
stage channel estimation in time domain is achieved through executing an inverse 
fourier transformation of received pilot symbols, executing a finite impulse re- 
sponse filtering of the inverse fourier transformation result, zero padding, and 
executing fourier transformation of the result of zero padding to yield tentative 
estimates (Z^) in relation to antennas in transmission antenna subsets (A^). 

30. Method according to claim 14 or 23, characterized in that the sec- 
ond stage channel estimation in time domain is achieved through executing an 
inverse fourier transformation of tentative estimates (Z^) in relation to anten- 
nas in transmission antenna subsets (A^), optionally executing a finite impulse 
response filtering of the inverse fourier transformation result, zero padding, and 
executing a fourier transformation of the result of zero padding to yield estimates 
of the channel reponse (Hffi. 

3 1 . Method according to claim 29 or 30, characterized in that it further 
comprises a step of interpolating the result of finite impulse filtering before fourier 
transformation thereof for synthesizing a number of tentative estimates in the time 
domain corresponding to the number of sub-carriers of the two dimensional data 
sequence. 

32. Method according to claim 29 or 30, characterized in that a ratio 
of the number of tentative estimates in the time domain to the number of received 
pilot symbols corresponds to pilot spacing in the first dimension of transmission. 

33. Method according to one of the claims 29 to 32, characterized in 
that the second stage channel estimation is excuted on the result of finite impulse 
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filtering. 

34. Method according to one of the claims 1 to 33, characterized in 
that it is applied to a celluar communication system with a frequency reuse factor 
of one such that base stations and related anntenna arrays form the plurality of 
transmit antennas and such that transmission antenna subsets CA Ml ) and related 
transmission antennas are defined in relation to this plurality of transmit antennas. 

3 5 . Channel estimator for estimating multiple input multiple output trans- 
mission channels in two dimensions, wherein transmission signals impinge from 
a plurality of (N T ) transmit antennas and carry a two dimensional data sequence 
with embedded pilot symbols, comprising: 

a memory storing a division of the plurality of transmit antennas into disjoint 
transmission antenna subsets G4 Ml ); 

at least one first estimator stage adapted to separate impinging pilot sequences 
in relation to transmission antenna subsets CA W ) and to perform a first stage chan- 
nel estimation to yield tentative estimates of a channel response in a first dimension 
of transmission; 

at least one second estimator stage adapted to seperate impinging pilot se- 
quences in relation to antennas in transmission antenna subsets (A^) and to per- 
form a second stage channel estimation for each antenna in each transmission an- 
tenna subset (y4 Ml ) to yield an estimation of the channel response. 

36. Channel estimator according to claim 35, characterized in that the 
first estimator stage is adapted to perform the first stage channel estimation using 
pilot sequences arranged as a two dimensional grid of pilot symbols, wherein pi- 
lot symbols used for first stage channel estimation depend on the first dimension 
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of transmission only and pilot symbols used for second stage channel estimation 
depend on a second dimension of transmission only. 

37. Channel estimator according to claim 36, characterized in that the 
first stage estimator stage is adapted to perform the first stage channel estimation 
using pilot sequences using pilot sequences expressed in a product form for achiev- 
ing seperability of pilot sequences in the first dimension of transmission and the 
second dimension of transmission. 

38. Channel estimator according to claim 37, characterized in that the 
first stage estimator stage is adapted to perform the first stage channel estimation 
using pilot sequences selected from orthogonal designs. 

39. Channel estimator according to claim 38, characterized in that the 
first stage estimator stage is adapted to perform the first stage channel estimation 
using Walsh pilot sequences or phase shifted pilot sequences. 

40. Channel estimator according to one of the claims 35 to 39, charac- 
terized in that the first estimator stage is adapted to perform the first stage channel 
estimation in frequency dimension. 

41. Channel estimator according to claim 40, characterized in that the 
first estimator stage is adapted to perform the first stage channel estimation using 
pilot sequences expressed in the product form according to 

= • X™ ' » = V2 + N T2 - (ft, - 1) 

Pi = {!,•■• ,N T i} 
ti 2 = {!,■■■ ,N T2 } 
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wherein 

I referes to a vector of data symbols in the two dimensional data sequence ar- 
ranged in time direction; 

% is a sub-carrier index in freqeuncy direction; 

represent pilot symbols for the first stage channel estimation in the first di- 
mension of transmission; 

X2? represent pilot symbols for the second stage channel estimation in a sec- 
ond dimension of transmission; and 

Nti is the number of transmission antenna subsets (A^), N T2 is the number of 
antennas in the transmission antenna subsets (y4. w ), and a product of the num- 
ber of transmission antenna subsets and the number of antennas in the transmis- 
sion antenna subsets equals the number Nt of the plurality of transmit antennas 
N T = N T i • N T2 . 

42. Channel estimator according claim 40 or 41, characterized in that 
the first stage estimator and/or the second stage estimator are adapted to achieve 
channel estimation for a pilot spacing having a value of one in a non-interpolating 
manner through yield of tentative estimates in relation to pilot symbol grid posi- 
tions. 

43. Method according to one of the claims 40 to 41, characterized in 
that the first stage estimator and/or the second stage estimator are adapted to 
achieve channel estimation for a pilot spacing having a value larger than one chan- 
nel estimation in an interpolating manner through yield of tentative estimates for 
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all data sequence grid positions in the dimension of estimation. 

44. Channel estimator according to one of the claims 41 to 43, charac- 
terized in that the first estimator stage is adapted to perform the first stage channel 
estimation using the product form of pilot symbols and a representation of received 
pilot symbols according to 

Y n = ]T X[f ]T X { 2 f B.f. + NT2 ■ (M1 - 1)} +% 
wherein 

2ppo A ^ ^ X^ 2 ^ H^ 2 + Nt2 ' ^ ~ 1 ' ) 

is a representation of a received pilot sequence to be estimated in relation to the 
transmission antenna subset ; 

H^j is the frequency response of transmit antenna /j; and 

Nl-- is noise superimposed during transmission. 

45. Channel estimator according to one of the claims 40 to 44, charac- 
terized in that the first estimator stage is adapted to perform the first stage channel 
estimation in frequency domain. 

46. Channel estimator according to one of the claims 40 to 44, charac- 
terized in that the first estimator stage is adapted to perform the first stage channel 
estimation in time domain. 

47. Channel estimator according to claim 45 or 46, characterized in that 
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second estimator stage is adapted to perform the second stage channel estimation 
in frequency domain. 

48. Channel estimator according to claim 45 or 46, characterized in that 
second estimator stage is adapted to perform the second stage channel estimation 
in time domain. 

49. Channel estimator according to one of the claims 35 to 39, charac- 
terized in that the first estimator stage is adapted to perform the first stage channel 
estimation in time dimension. 

50. Channel estimator according to claim 49, characterized in that the 
first estimator stage is adapted to perform the first stage channel estimation using 
pilot sequences expressed in the product form according to 

X% = X$>-X%>\ f x = ^ + N T2 -(i, 1 -l), ^ = {1,--- ,N T1 j 

to = {l,--- ,N T2 } 

wherein 

I referes to a vector of data symbols in the two dimensional data sequence ar- 
raged in time direction; 

i is a sub-carrier index in freqeuncy direction; 

represent pilot symbols for the second stage channel estimation in the first 
dimension of transmission; 

X£ 2) represent pilot symbols for the first stage channel estimation in a second 
dimension of transmission; and 

Nt\ is the number of transmission antenna subsets (*4 Ml ), iV T2 is the number of 
52 
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antennas in the transmission antenna subsets (A^), and a product of the num- 
ber of transmission antenna subsets and the number of antennas in the transmis- 
sion antenna subsets equals the number N T of the plurality of transmit antennas 
N T = N T i • N T2 . 

51. Method according to one of the claims 49 or 50, characterized in 
that for a pilot spacing having a value of one first stage channel estimation and/or 
second stage channel estimation is achieved in a non-interpolating manner through 
yield of tentative estimates in relation to pilot symbol grid positions. 

52. Method according to one of the claims 49 or 50, characterized in 
that for a pilot spacing having a value larger than one first stage channel estima- 
tion and/or second stage channel estimation is achieved in an interpolating manner 
through yield of tentative estimates for all data sequence grid positions in the di- 
mension of estimation. 

53. Channel estimator according to one of the claims 50 to 52, charac- 
terized in that the first estimator stage is adapted to perform the first stage channel 
estimation using the product form of pilot symbols and a representation of received 
pilot symbols according to 




wherein 

N T2 

1*2 = 1 

is a representation of a received pilot sequence to be estimated in relation to trans- 
mission antenna subset A^ ; 
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Hfj is the frequency response of transmit antenna fi; and 
Nfj is noise superimposed during transmission. 

54. Channel estimator according to one of the claims 49 to 53, charac- 
terized in that the first estimator stage is adapted to perform the first stage channel 
estimation in frequency domain. 

55. Channel estimator according to one of the claims 49 to 53, charac- 
terized in that the first estimator stage is adapted to perform the first stage channel 
estimation in time domain. 

56. Channel estimator according to claim 54 or 55, characterized in 
that the second estimator stage is adapted to perform the second stage channel 
estimation in frequency domain. 

57. Channel estimator according to claim 54 or 55, characterized in 
that the second estimator stage is adapted to perform the second stage channel 
estimation in time domain. 

58. Channel estimator according to claim 45 or 54, characterized in that 
the first stage estimator comprises a first finite impulse response filter adapted to 
perform the first stage channel estimation in frequency domain is through finite 
impulse response filtering of received pilot symbols to yield tentative estimates 
(Z^) of received pilot sequences (Zff) in relation to transmission antenna sub- 
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sets (Am). 

59. Channel estimator according to claim 47 or 56, characterized in 
that the second stage estimator comprises a second finite impulse response fil- 
ter adapted to perform the second stage channel estimation in frequency domain 
through finite impulse response filtering of tentative estimates (Zff) in relation to 
antennas in transmission antenna subsets (A^) to yield estimates of the channel 
reponse (Hfl). 

60. Channel estimator according to claim 58 or 59, characterized in 
that the second stage estimator comprises a buffer for executing the second stage 
channel estimation in frequency domain through finite impulse response filtering 
with induced time delay for smoothing type filtering. 

6 1 . Channel estimator according to one of the claims 58 to 60, character- 
ized in that the first finite impulse response filter and/or the second finite impulse 
response filter is a Wiener type finite impulse response filter. 

62. Channel estimator according to claim 46 or 55, characterized in that 
the first stage estimator comprises: 

an first inverse fourier transformation unit adapted to execute an inverse fourier 
transformation of received pilot symbols; 

a third finite impulse response filter adapted to execute an optional finite im- 
pulse filtering of the inverse fourier transformation result; 

a zero padding unit; and 
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a first fourier transformation unit adapted to execute a fourier transformation of 
the output of the zero padding unit to yield tentative estimates (Zfl ] ) in relation to 
antennas in transmission antenna subsets (A^) as first stage channel estimation in 
time domain. 

63 . Channel estimator according to claim 48 or 57, characterized in that 
the second stage channel estimator comprises: 

a second inverse fourier transformation unit adapted to execute an inverse fourier 
transformation of tentative estimates (Zff) in relation to antennas in transmission 
antenna subsets (A^); 

a fourth finite impulse response filter adapted to execute an optional finite im- 
pulse filtering of the inverse fourier transformation result; 

a zero padding unit; 

a second fourier transformation unit adapted to execute a fourier transformation 
of the output of the zero padding unit to yield estimates of the channel reponse 
{Hfl) as second stage channel estimation in time domain. 

64. Channel estimator according to claim 62, characterized in that it 
further comprises an interpolation filter adapted to interpolate the result of finite 
impulse response filtering before fourier transformation thereof for synthesizing a 
number of tentative estimates in the time domain corresponding to the number of 
sub-carriers of the two dimensional data sequence. 

65. Channel estimator according to claim 64, characterized in that a 
ratio of the number of tentative estimates in the time domain to the number of 
received pilot symbols corresponds to pilot spacing in the first dimension of trans- 
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mission . 

66. Channel estimator according to one of the claims 62 to 65, charac- 
terized in that the second stage estimator is connected to the third finite impulse 
response filter and adapted to process the result of related finite impulse filtering 
for determination of the second stage channel estimation. 

67. Channel estimator according to one of the claims 35 to 66, charac- 
terized in that it is operable in a celluar communication system with a frequency 
reuse factor of one where base stations and related anntenna arrays form the plu- 
rality of transmit antennas and such that transmission antenna subsets (A^) and 
related transmission antennas are defined in relation to this plurality of transmit 
antennas. 



68. Computer program product directly loadable into the internal mem- 
ory of a channel estimator for estimating multiple input multiple output transmissio 
channels in two dimensions comprising software code portions for performing the 
steps of one of the claims 1 to 34 when the product is run on a processor of the 
channel estimator for estimating multiple input multiple output transmissio chan- 
nels in two dimensions. 
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